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Preface 

These three reports describe selected applications that the Texas Instruments Digital Signal Processing 
(DSP) Group, in conjunction with TI research groups, has developed in the emerging field of fuzzy logic. 

The first paper gives an overview of the fuzzy logic theory and suggests methods of implementation that 
are illustrated in the next two papers. 

The second paper describes the implementation of fuzzy logic as a serial algorithm on a TMS320 DSP sys- 
tem. The algorithm is based on parts of the theory described in the third paper. 

The third paper presents a theory for a dedicated silicon-level fuzzy logic implementation. 
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What is Fuzzy Logic? 
An Overview of the 
Latest Control Methodology 



Timothy A. Adcock 
Digital Signal Processing — Semiconductor Group 
Texas Instruments Incorporated 



Introduction 



The name "Fuzzy Logic" seems to imply an imprecise methodology that is useful only when accuracy is 
not necessary or important. That is what many people assume when they first hear about fuzzy logic — and 
understandably so. In a world increasingly manipulated by computers with their absolute "1" or "0" and 
"on" or "off' concepts, a term like fuzzy logic suggests inaccuracy or imprecision. Even Webster's dictio- 
nary defines "fuzzy" as: 

fuzz«y (-e) adj. 2. not clear, distinct, or precise; blurred 

This is not true of fuzzy logic. Fuzzy logic can address complex control problems, such as robotic arm 
movement, chemical or manufacturing process control, antiskid braking systems, or automobile transmis- 
sion control with more precision and accuracy, in many cases, than traditional control techniques have. 

Fuzzy logic was invented and named by Lotfi Zadeh, a professor at the University of California at Berkeley. 
Fuzzy logic is a methodology for expressing operational laws of a system in linguistic terms instead of 
mathematical equations. Many systems are too complex to model accurately, even with complex mathe- 
matical equations, but fuzzy logic's linguistic terms provide a useful method for defining the operational 
characteristics of such a system. These linguistic terms are most often expressed in the form of logical 
implications, such as If - Then rules: 

Ifairjemp is WARM, then set fan_speed to MEDIUM. 

The terms WARM and MEDIUM are actually sets that define ranges of values known as membership func- 
tions. By choosing a range of values instead of a single discrete value to define the input variable 
"air_temp", you can control the output variable "fan_speed" more precisely. Fuzzy logic controllers can 
often improve the performance of a control system by reducing the chance of wild functions in the output 
that may be caused by variations in the measured input variables. 

The Traditional Approach 

To illustrate the difference between fuzzy logic and the traditional approach, here is a control problem. 
First, consider how the traditional — often called "crisp" — controller would handle it: 

Ifairjemp is £ 70° Fahrenheit, then set fan_speed to "1000 rpm". 
Ifairjemp is < 70° Fahrenheit, then set fan_speed to "100 rpm". 

A nonfuzzy, or "crisp," controller relies on a discrete valued decision point. For this type of system, the 
input must reach an exact value before the control system reacts in a certain way. Even small variances in 
this input value may cause the output to react drastically differently. For instance, if the temperature is 70° 
or above, the first rule will set the fan_speed to 1000 rpm. If the temperature is below 70°, the second rule 
will set the fan_speed much lower, to 100 rpm. Figure 1 shows a diagram of this crisp valued controller. 
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Figure 1. Crisp Controller 
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What would happen if the temperature were 69.5 ° ? Or, more importantly, what would happen to the control 
system if the temperature were transitioning from below 70° to above 70°? The temperature might even 
fluctuate back and forth slightly above or below 70° (e.g., 69.0° to 71.0°). This would cause the control 
system to alter the fan speed wildly for changes in the input variable air_temp, although the temperature 
change may not be significant. 

These transition points are difficult for "crisp" control systems to handle, but they are exactly where "fuzzy 
logic" excels. 

Fuzzy Control 

Fuzzy logic is implemented in three phases (see Figure 2): 

1. Fuzzification (crisp input to fuzzy set mapping). 

2. Inference (fuzzy rule generation). 

3. Defuzzification (fuzzy to crisp output transformation). 

Figure 2. Fuzzy Logic Phases 
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Fuzzification 

In the first fuzzy logic phase — fuzzification — actual measured input values are mapped into fuzzy mem- 
bership functions. As an example, a climate-control system has been developed with fuzzy logic. 

To create a climate control system, we first developed membership functions for the input variable 
"airjemp". These membership functions are defined by both a range of values and a degree of member- 
ship. In fuzzy logic, it is important to distinguish not only which membership functions a variable belongs 
to, but also the relative degree to which it is a member. This gives the variable a "weighted" membership 
in a membership function. A variable can have a weighted membership in several membership functions 
at the same time. The membership functions for "air_temp" are shown in Figure 3. 
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Figure 3. Air_Temp Input Variable Membership Functions 




As shown in Figure 3, fuzzy membership functions span a range of values and can actually overlap. Three 
sets of membership values are defined above for the variable "air_temp". They are COOL, WARM, and 
HOT. The degree of membership is found by finding the intersection point of a distinct input value on the 
horizontal axis with the line defining one or more fuzzy membership functions. This intersection point is 
assigned a corresponding value on the vertical axis to define the relative membership in a set for an actual 
measured input value. Notice that when "air_temp" is at a particular value, it may be contained in one or 
more fuzzy sets. For instance, at 70° , "air_temp" is a member of the function HOT with a relative member- 
ship of 0.17. It is also a member of the function WARM with a relative membership of 0.37. Unlike a crisp 
system in which a value either is or is not a member of a function, a fuzzy logic system can take action based 
not only on membership in a fuzzy set, but also on the degree to which a variable is included in a member- 
ship function. In this case, because "airjemp" at 70° is more WARM (0.37) than it is HOT (0.17), the 
controller will take that into account when defining what output action to take. 



Inference Rule Definition 



Once membership functions have been defined for input and output variables, a control rule base can be 
developed to relate the output actions of the controller to the observed inputs. This phase is known as the 
inference, or rule definition portion, of fuzzy logic. Any number of rules can be created to define the 
actions of the fuzzy controller. Some examples are shown below. 



Fuzzy Logic Rule Definition 

If air_temp is COOL, then set fan_speed to SLOW. 

If airjemp is HOT, then set fan_speed to FAST. 

If airjemp is WARM, then set fan_speed to MEDIUM. 
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These If-Then rules can relate multiple input and output variables. Because the rules are based on word 
descriptions instead of mathematical definitions, any relationship that can be described with linguistic 
terms can typically be defined by a fuzzy logic controller. This means that even nonlinear systems can be 
described and easily controlled with a fuzzy logic controller. In addition, since variables have weighted 
memberships — in particular membership functions — the rules that are composed of these variables are 
weighted as well. This means that different rules have different impacts on the controller, according to the 
measured input variable. For a multiple-input/multiple-output system with many defining rules, a wild 
fluctuation in any single input will be tempered by these rule weightings. Because of this, fuzzy logic sys- 
tems are very robust and often allow many rules to be removed or altered without significantly impacting 
the controller. 

Defuzzification 

After the fuzzy logic controller evaluates inputs and applies them to the rule base, it must generate a usable 
output to the system it is controlling. This may mean setting a voltage or current to a particular value to 
control the speed of a fan in the example above, or it may mean defining the optimal speed of a robotic arm 
as it nears its target. The fuzzy logic controller must convert its internal fuzzy output variables into crisp 
values that can actually be used by the controlled system. You can perform this portion of the fuzzy control 
algorithm, known as defuzzification, in several ways. Two of the most common methods are: 

• maximum defuzzification method (page 6). 

• centroid calculation defuzzification method (page 7). 

Remember from fuzzification that in mapping input variables to membership functions, a particular mea- 
sured value of the input variable determined the relative membership of that input variable in an input mem- 
bership function. To determine the mapping of output variables to their corresponding output membership 
functions, the weighted input membership function and corresponding rule base determine the relative 
membership in the output function. Whatever relative membership was given to the input variable will also 
be given to the output variable, as assigned by its corresponding rule. For air_temp = 70°, the output vari- 
ables are assigned a value that corresponds to the input value shown in Table 1. 

Table 1. Fan_Speed (Membership Function Relative Membership) 



Input Variable 


Defining Rules 


Output Variable 


airjemp (WARM) = 0.37 


If airjemp = WARM, 

then set fan_speed to MEDIUM 


fan_speed (MED) = 0.37 


airjemp (HOT) = 0.17 


If airjemp = HOT, 

then set fan_speed to FAST 


fan_speed (FAST) = 0.17 



The output variable fan_speed is given the same relative mapping as the input variable airjemp that is 
defined by a particular rule. 

Figure 4 illustrates the output variable membership functions. It this case, a distinct value on the horizontal 
axis is defined by the relative membership on the vertical axis. To create the actual crisp output value for 
the controller system output, membership functions are used with the output variable. In Table 1, the input 
value airjemp = 70 resulted in two weightings for fan_speed: 

• Fan_speed = 0.37 was assigned to the output membership function MEDIUM. 

• Fan_speed = 0.17 was assigned to the output membership function FAST. 
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As shown in Figure 4, the actual output value is determined by beginning at the weighting factor on the 
vertical axis and moving horizontally until an intersection point is reached on the lines defining its 
associated membership function. This intersection point is then transposed to the horizontal axis to deter- 
mine the crisp output value. 




Fan Speed (rpm) 



Maximum Defuzzificatlon Method 

One method of defuzzification is known as the maximum method. In this method, if more than one rule 
is active, the maximum relative membership is used to determine the output value. In the above example, 
making air_temp = 70° created two possible values for fan_speed: 

fanjspeed = 0.37 
fan_speed = 0.17 

The maximum defuzzification method provides a single output by choosing the active rule with the greatest 
relative membership value in the output membership function. In the preceding example, the following rule 
is chosen because it has the highest membership value for fan_speed. 

Ifairjemp = WARM, then set fanjspeed to MEDIUM. Fanjspeed (MED) = 0.37 

The value 0.37 on the vertical axis intersects the membership function MEDIUM at two points — one on 
the positive slope (at 375 rpm) and one on the negative slope (at 710 rpm) of the function. The two points 
represent two possible solutions that must be resolved. 
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Centroid Calculation Defuzzification Method 



Another method for calculating the output value is the centroid method. In this method, a weighted average 
of all the active rules determines an output by summing all of the applicable output variables over their 
relative membership values. Although this method is more computationally intensive, it creates a distinct 
output value based on the relative memberships of all of the active rules that apply (see Figure 5). This 
method eliminates the problem of multiple solutions observed with the maximum method. A processor 
architecture with a hardware multiply-accumulate feature like that of the TMS320 DSP family excels at 
this method. 

Figure 5. Fan_Speed Output Centroid Calculation 




Centroid Point 
(635 rpm) 

Fan Speed (rpm) 



Summary 

By using fuzzy logic, you can simplify complex control problems that once required a high-powered micro- 
processor to execute in real time; you can now execute them on a low-cost Texas Instruments TMS320 DSP 
or TMS370 microprocessor. The following application note shows the benefits of controlling a simple DC 
motor with fuzzy logic using a TMS320C14 digital signal processor. 



Implementation of 
Fuzzy Logic Servo Motor Control 

on a 
Programmable 
Texas Instruments TMS320C14 DSP 



Mathew George, Jr. 
Digital Signal Processing — Semiconductor Group 
Texas Instruments Incorporated 



Abstract 



This paper describes the implementation of a fuzzy logic compensator on a Texas Instruments TMS320C14 
DSP-based servo motor control development system. The system contains a real motor that is controlled 
by the programmable DSP. An on-chip debugger and servo motor program allowed both simple code modi- 
fication and interactive control of the motor. A fuzzy logic algorithm was directly substituted for the origi- 
nal PID algorithm; this resulted in comparable motor response and algorithm performance. This imple- 
mentation proves the feasibility of real-time fuzzy logic-based servo motor control on a real system. 



Introduction 

Fuzzy logic is relatively new theory. Most of the readily available hands-on fuzzy logic system examples 
have been software simulations or bulky real systems. However, a TI commercial microprocessor (the 
TMS320C14) can serve as a simple, real-time, real-system platform for applying and investigating fuzzy 
logic. This facilitates both the understanding and implementation of fuzzy logic as a real-time program- 
mable solution for the general engineering public. 

Servo motor control is a viable and useful implementation. A programmable PID motor control board uses 
a Texas Instruments TMS320C14 chip and has an actual motor whose performance can be observed. Faster 
and newer parts are available, but the 'C14 is optimized for motor control with such features as on-board 
pulse-width modulation (PWM) generation capabilities. The board also has on-chip debugger code and 
interactive PID control code. The PID compensator code (written in TMS320 assembly language, which 
is upward-compatible with code executed by such newer fixed-point TI DSPs as the TMS320C25 or 
TMS320C50) uses position and velocity of the motor for inputs and motor input current as the output. The 
control code will allow the operator to interactively change such values as servo position and velocity and 
to monitor position error. Modifying the code for fuzzy logic required replacing the PID compensator sec- 
tion of the code with a fuzzy logic compensator. 

The membership function for the compensator defines the error between present motor position and the 
desired (command) position of the controller. Five linguistics variables characterize the function: negative 
medium, negative small, zero, positive small, and positive medium. The function is represented as overlap- 
ping isosceles right triangles for ease of fuzzification. Various rules for an inverted pendulum control sys- 
tem (ball and stick) are explained in [1]. These same eleven rules were used for servo motor control. 

The algorithm's three sections are implemented as a series of software loops: fuzzification (i.e., input evalu- 
ation), fuzzy inference (rule contribution that uses a table look-up), and defuzzification (using center-of- 
gravity method). It is based on an algorithm developed for [4]. Each section uses various arrays that are 
modified in that section. The sizes of these loops are directly proportional to the number of inputs, outputs, 
and rules used in the system. 



Servo Motor System (Power-14) 

The heart of the Power-14 board is a Texas Instruments TMS320P14 chip (one-time programmable 
TMS320C14). (See Figure 1). 
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Figure 1. Power 14 System 




You communicate to the board through an RS-232 serial port connection by using standard terminal or ter- 
minal emulation software (such as Procomm). The 'P14 peripherals are optimized for control applications. 
An event manager can be operated in a PWM (pulse-width modulation) mode that is ideal for motor control. 
Monitor code is loaded into 'P14 external memory and run, which provides a command line debugger. The 
debugger has all standard debugger functions, such as memory dumps and modification, stepping through 
code, breakpoints, etc. The monitor can be used to load and run the servo motor program with PID compen- 
sator. The program is interactive and lets you control the motor from the keyboard. Position, velocity, and 
the PID values can be set. Data acquisition functions allow an ASCII text input stimulus table to be loaded 
onto the Power-14, an acquisition run to be executed, and the resulting ASCII output table to be written 
to a PC file for graphing. The rest of the board contains support circuitry: amplifiers for the motor and a 
serial port interface. An encoder on the motor is used as a position sensor for a compensator input. The 
velocity is found by executing a back-difference. Note that there is no separate velocity sensor. 



PID Implementation 

The source code for this system implements the PID compensator in one file (See Appendix A). The algo- 
rithm is a direct implementation of the PID equation. The Proportional, Integral, and Differential variables 
are derived from the motor encoder sensor detecting position. On each cycle, the present position is taken 
from the encoder and stored in Position. The error is found by subtracting Position from DesiredPosition 
and storing it in ErrNow. Thus, ErrNow is the position input for the proportional section of the PID . A back- 
difference is then taken with ErrNow and ErrLast (the error in the previous cycle) to give ErrDiff . ErrDiff 
is an approximation of the velocity and is therefore the second input (Differential) of the PID. The Integral 
is found by adding the ErrNow value and storing it in Kintegrator. The following equation then holds the 
compensator output: 

PWM = 9* P* ErrNow + 9 * I * Kintegrator + 9* D* ErrDiff a) 

The value is scaled for the PWM mode and stored in NewServo. Thus, for the compensator output, the 
actual motor current input value is converted to a PWM value for implementation. The PWM output is then 
sent to power amplifiers that finally drive the motor. The PWM frequency can be controlled from the PID 
program. 
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Fuzzy Logic Theory for Servo Motor Control 



The membership function for this system is simple. Five fuzzy logic ranges (linguistic variables) were cho- 
sen to remain consistent with [1] and used for both inputs. Sets of values are represented as overlapping 
isosceles right triangles (Figure 2). 

Figure 2. Membership Function 
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Key 

PM=positive medium 
PS=positive small 
Z=zero 

NS=negative small 
NM=negative medium 

Two input variables, position (Theta) and velocity (dTheta) of the motor, are used in this fuzzy logic system. 
One output variable, motor current, is used, which will be proportional to the PWM output that is actually 
written. The rules for the compensator as mentioned were taken from [1]. Thus Theta, dTheta, and motor 
current operate according to the following "if a and b, then c" rules, as shown in Table 1. 

Table 1. List of Rules 



If Theta = 


And dTheta = 


Then Motor Current = 


Z 


Z 


Z 


PS 
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NS 


PM 
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NM 


NS 
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PS 


NM 
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These rules can then be indexed according to the scale shown in Figure 2. The mapping seen in Table 2 
will be used in TMS320 programming. 

Table 2. Indexed List of Rules 



If Theta = 


And dTheta = 


Then Motor Current = 











1 





-1 


2 





-2 


-1 





1 


-2 





2 





-1 


1 





-2 


2 
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-1 
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-2 


1 


-1 





-1 


1 






The defuzzification is done by the center-of-gravity method. 

Fuzzy Logic Implementation 

Arrays are used in the fuzzy logic calculations and modified in the various loops that implement the com- 
pensator. Appendix B lists the TMS320C14 code. Figure 3 shows the three sections of the compensator: 
fuzzification, fuzzy inference, and defuzzification. Notation for the arrays follows C language standard, 
with subscripts from to n-1. The 'C14 algorithm is based on an algorithm developed for [4]. The figure 
key summarizes the values of the system that will be used in the examples in this section. 
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Figure 3. List of Arrays 
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X[ninputs][nmembs] 


Fuzzy 
Inference 


for i=0 to nrules-1 
for j=0 to ninputs-1 
for k=0 to noutputs-1 








U[noutputs][nmembs] 


Defuzzification 


for i=0 to nmembs-1 
for i=0 to nmembs-1 




1 


u [outputs] 





Fuzzy Logic Implementation as a Series of Loops 



Key 
nmembs=5 
ninputs=2 
nrules=11 
noutputs=1 

i, j, k = loops within loops 

Note which arrays are modified in each section and their size boundaries. They are discussed later in more 
detail. The two inputs seen in x[ninputs] are position (found from the encoder) and velocity (found from 
an approximation of the derivative by taking the back-difference of the position). These are the Theta and 
dTheta variables described previously. As the compensator code begins, the position and velocity inputs 
(ErrNow and ErrDiff, respectively) are copied to the array X[ninputs]. The position is mapped so that one 
rotation of the motor ranges from -255 to +255 (See Figure 4 a). This relationship is mapped onto the x-axis 
of the membership function and thus fuzzifies the position of the motor (Figure 4 b). Note that this figure 
is not drawn to scale. 



14 



Figure 4. Motor Shaft/Membership Function Mapping 

a. Motor Range 
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In the fuzzification loop, the degree of membership of each input relative to the input membership function 
is evaluated and written in the array X[ninputs][nmembs]. The value for a particular linguistic variable is 
the y value of the triangle for a particular value of x (See Figure 5). 

Figure 5. Analysis for One Linguistic Variable 
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The y value is found by using the simple algebraic equation for a line (y=mx+b). This equation may be 
geometrically reduced to one of the two following equations, depending on which side of the triangle the 
value of x lies: 

if {x : xO < x < xl} then y = (x - xO) / (xl - xO) else (2) 
if{x:xl<x< x2} then y = (x2 - x) / (x2 - xl) 

The division needed is costly on most microprocessors, usually requiring at least a number of cycles equiv- 
alent to the number of bits of the number being divided. But if the slope m is made equal to 1 by causing 
the elements of the membership functions to be isosceles right triangles, the equation can be reduced to 
y=x-xO. This translates into a simple one-cycle subtraction. Note that forcing the elements of the member- 
ship functions to be isosceles right triangles also forces the peak of the membership function to no longer 
be 1. Rather the peak value=xl-x0=x2-xl=85 (as seen in Example 1). This action also eliminates the need 
for using a Q format [6] to represent the fractional values from Equation 2 if the triangle were not isosceles. 



Example 1 demonstrates this fuzzification calculation. If the position input value were 100(i.e.,x[0]=100), 
it would have nonzero degrees of membership in the PS and PM linguistic variables. You can also see this 
in Figure 5. To calculate the actual degree of membership value, the value for x is plugged into Equation 
2 for both PS and PM boundaries. Thus, in PS the contribution is 70, while in PM it is 15. The rest of the 
linguistic variables are because there is no contribution. (Of course, in software, all linguistic variables 
must be evaluated.) For this example, X[0][nmembs]=[0, 0, 0, 70, 15]. 

Example 1. Fuzzification Example 




-255 



x=100, thus 
y= 70 and y =15 
PS PM 



The next step involves fuzzy inference. In this loop, the maximum and minimum functions, as explained 
in [2], are implemented. In the actual code, the array indexing of the membership values is made nonnega- 
tive by adding a bias of three. Therefore, instead of NM to PM being indexed from -2 to +2, as seen in 
Figure 4, they are indexed from to 4. Table 3 shows how the 11 rules are indexed and reindexed in the 
array RULE_TABLE[rule] [input+output] . 
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Table 3. Original and Reindexed Table of Rules 



Original Index 


Reindexed 


x(0) 


x(l) 


u(0) 


x(0) 




u(0) 











2 


2 


2 


1 





-1 


3 


2 


1 


2 





-2 


4 


2 





-1 





1 


1 


2 


3 


-2 





2 





2 


4 





-1 


1 


2 


1 


3 





-2 


2 


2 





4 





1 


-1 


2 


3 


1 





1 


-2 


2 


4 





1 


-1 





3 


1 


2 


-1 


1 





1 


3 


2 



Some explanation is required for this decoding. The indexed rules match the explicit rules as described in 
Table 1. Also, the values given by accessing the RULE_TABLE are limited to the values indexed by 
nmembs. This characteristic is heavily used in the index manipulation and allows nmemb to be inter- 
changed with RULE_TABLE[rule][input+output] in the appropriate parts of the algorithm. 

To find the minimum value of X[ninputs] [nmembs] decoded from the rule table inputs and stored in minZ 
(which is initialized to the maximum y value — in this case, 85), the equation is: 

minZ = min (X[input] [RULEJTABLE [rule] [input]]) p) 

(Since noutputs=l, the loop is simplified, and minZ does not need to be the general case array 
minZ[nrules]). Note that only the first two columns (the input columns) of RULEJTABLE are used in this 
part of the fuzzy inference section. 

Then, for each rule, the max is taken of the output value U[nmembs], which is initialized to 0. The general 
case U[output][ nmembs] is simplified because only one output is decoded from the rule table outputs and 
minZ. This equation can be summarized as: 



U[nmembs] = U(RULE_TABLE[rule] [ninput + output]) 

= max (U[RULE_TABLE[rule] [ninput + output]], mini) 




Note that in this section only the last column (the output column) is used. The following equation summa- 
rizes the minimum and maximum functions that are executed for each rule to result in the array U[nmembs] 
by plugging Equation 3 into Equation 4: 

U[nmembs] = max {U[RULE_TABLE[rule] [ninput + output]] , (J) 
min (X[input] [RULE_TABLE[rule] [input]])} 
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The following example illustrates the fuzzy inference section. One cycle of the loop for rule=0 is shown. 
For the input array: 



X[ninputs] [nmembs] = I 70 15 | 
= | 60 25 | 



(6) 



and for rule=0, input=0, and input=l, plug into Equation 3: 

minZ = min | (X[0] [2]) | = min | 70 | =25 
| (X[l] [2]) | | 25 | 



(7) 



then for the max rule=0, input=0, and input=l, plug into Equation 4: 

U[nmembs] = U[2] = max | U{2] \ = max | | = 25 

| mini | I 25 | 



(8) 



This process continues for the list of 11 rules so that the array U[nmembs] contains the maximum values 
of the min — i.e., the contribution of each rule to the inference. 

The final step involves defuzzifying the U[noutputs] [nmembs] array. Since noutputs=l, U simplifies to 
U[nmembs]. The U[nmembs] array now has five values in it for its corresponding five positions. The cen- 
ter-of-gravity calculation is done with two loops. The first finds the numerator by using the multiplier to 
weight Ufnmembs] by its position. The second loop finds the denominator by summing the position. The 
inevitable 16-bit divide loop then finds the output value u[output]. This u[output] represents the motor cur- 
rent mentioned in Fuzzy Logic Theory for Servo Motor Control (page 12). The divide operation is done on 
the 'C14 by a 16-cycle loop. This value is then scaled for the output and written to the NewServo memory 
location that sends it to the PWM generator. 



The fuzzy logic code with two inputs, eleven rules, one output, and five linguistic variables requires about 
2000 instruction cycles to execute on the TMS320C14. This means a 400-us period (2.5-kHz frequency) 
because one instruction cycle is 200 nanoseconds on the 'C14. The update period of the motor is 3.8 ms 
(263-Hz frequency), so the fuzzy code is obviously adequate (for comparison, the PID code execution 
required nine (is (111-kHz frequency), for the update period). Figure 6 shows the PID performance for a 
step with the classical control overshoot and settling times. 



As an example of defuzzification, 
if U[nmemb] = [0 15 70 35 0], 
then the output is: 

U[output] = 0*(-170) + 15*(-«5) + 70*(0) + 35*(85) + 0*(170) 



= 14.17 



These loops thus evaluate the control input necessary for the servo motor. 



Results 
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Figure 6. PID Controller 
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The fuzzy logic curve (Figure 7) is smoother than the PID controller curve but doesn't go to zero 
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Figure 7. Fuzzy Controller 
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The update period is adequate for applications such as servo motor, robotics, motion control, and automo- 
tive control. Better performance may be needed for such applications as hard-disk drives. Later-generation 
digital signal processors, such as the TMS320C5x, operate at up to 25 nanoseconds with much more effi- 
cient instruction code. 



Conclusion 

The final fuzzy logic system behaved favorably when compared to the conventional PID system. The sys- 
tem proved the feasibility of implementing a real-time fuzzy logic servo motor control. Proof of the fuzzy 
control was shown by positioning the motor spindle with an error outside the membership function, thus 
causing the control to desist. Further development could include a graphics display and the ability to vary 
the rules. The TMS320C14 board fits in a 12 x 8 x 6-inch suitcase conveniently and requires only an AC 
power supply and an RS-232 keyboard connection. This product easily demonstrates real-system fuzzy 
logic control on a microprocessor. 

Fuzzy logic has great potential as a programmable solution for general engineering. For applications where 
performance is a priority, a hard-wired silicon solution (which may even configure as a microprocessor 
peripheral) based on [4] is being developed. The programmable solution may assist in the transition to this 
hard-wired option, depending on the software/hardware tradeoffs. 
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APPENDIX A 
PID Code 



title "Servo compensator SRevision: 3.4 $" 

*••••••••••••••••*•***•**•*****•**••*•***•*•**••*************••••*••******* 

SHeader:: C:/src/cl4/ps/vcs/comp.s_v 3.4 01 Oct 1991 17:24:18 "$ 

$config$='7T8 /K! /L;*.ref.def /R:*- /B80" 
! config! ="/Mcomp.s" 

!NAME! 

comp.s 

1PATHS! 

modules 

!0! 

DESCRIPTION 

Servo compensator - uses a PID algorithm 

PRINCIPLE AUTHORS: 
Dave Sewhuk 

CREATION DATE: 

December 23, 1990 22:54:08 

COPYRIGHT NOTICE: 

(QCopyright 1990 Teknic Inc. All rights reserved. 

lend! 

*************************************************************************** 

********»♦*»»»***■,***********»*****♦*****»»**»******+»,*»*»****+*«*****»**+ 

HEADERS UTILIZED 

.include "macrodef.inc" 
.include "cl4io.inc" 

********************************************************* ****************** 

*********************************************************************** 

INAME! 

comp.s 



IPATHS! 
!0! 



Exported Variables 

.def PwmChannel,PwmPeriod 
.def NewServo,ErrDiff,Kintegrator 
.def ErrNow,ErrLast,DesiredPosition 



.bss PwmChannel.l 

.bss PwmPeriod,l 

.bss NewServo.l 

.bss ErrDiff,l 

.bss Kintegrator,l 

.bss ErrNow,l 

.bss ErrLast.l 

.bss DesiredPosition,! 



; Current PWM channel 

; Current period of PWM channels 

; New servo value 

; Error difference 

; Error integrator 

; Current error function 

; Last error function result 

; Servo to this position 



!end! 



*************************************************************************** 



1NAME! 

comp.s 

IPATHS! 

Local Variables 

!0! 
lend! 



•I************************************************************************* 
*#**************************************^ 



.text 

.def servo_ISR 

; 1NAME! 
servo_ISR: 

; IPATHS! 



This is the servo compensator. The distribution algorithm is 
the PID format. The compensator performs the following function 
PWM = 9*P*error + SUM(9*ek)*I + 9*(err diff)*D 
with the PWM output clipped to timer hardware limits. 



No changes to foreground environment. 



STACK LEVELS: 


EXTERNAL DATA REFERENCES 



.ref ONE.CONSTFFFF ; Useful constants 



functions\all 

functions\interrupt 

modules\utils.s 



!0! 

DESCRIPTION: 



RETURNS: 



.ref Position 



; Position from feedback ISR 
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.ref Pgain.Igain.Dgain ; Servo constants from user interface 
.ref ISR_TMP ; Scratch 



♦♦a****************************************************** 



!skip start! 



; Saturate math 



SOVM 

i 

LAC ErrNow 
SACL 



ZALH 
; Calc new error 

SUBH 
SACH 



En-Last 

DesiredPosition 

Position 
ErrNow 



Run PID algorithm on new sample 



LAC ErrNow 
SUB ErrLast 
SACL 

ZALH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
SACH 



ZALH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
SACH 

; Saved scaled up error 

ZALH 
; Update integrator 

ADDH 
SACH 



ErrDiff 

ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 

ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ISR TMP 



Kintegrator 

ErrNow 
Kintegrator 



; Save old error 



; Update error difference 



Scale up value a bit/saturation 

*2 

*3 

*4 

*5 

*6 

*7 

*8 

*9 



; Scale up/saturation 

;*2 
;*3 
;*4 
;*5 
;*6 
;*7 
;*8 

;*9 



LT Pgain 
MPY ISRTMP 
PAC 

LT Igain 
MPY Kintegrator 

LTA Dgain 
MPY ErrDiff 
APAC 
SACH 



NewServo 



; Do proportional part 

; Do integral part 
; Do differential part 
; Save result 



Check the answer so that it fits within the limits of the timer's 
period: [O..PwmPeriod*4) 



pwmisr20: 



LAC PwmPeriod,l 
SUB ONE.3 
SUB NewServo 
BGZ pwmisr20 
LAC PwmPeriod,l 
SUB ONE.3 

SACL NewServo 
B pwmisrSet 

LAC NewServo 

ADD PwmPeriod,l 

BGEZ pwmisrSet 

ZAC 

SUB PwmPeriod,l 

SACL NewServo 



Set Output PWM to new value 
pwmisrSet: 

LAC PwmPeriod,l 
ADD NewServo 
SACL NewServo 
.if ChipVlRl 
.ref CONST3 
SUB CONST3 
BLEZ pwmisrlO 
; No fixes needed 

LAC NewServo 
AND CONST3 
SUB CONST3 
BNZ pwmisrlO 
LACK 4 
ADD NewServo 



Check positive: 2*Period-4 
- value 

Oops, pin at biggest 

; Negative OK? 
; Add largest offset 

; Make largest negative 



; Get center period value 
; Add newly computed value 

; PWM bug workaround 



; Lower 2 bits set? 
; OK No fixes needed 
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SACL 



NewServo 



pwmisrlO: 



.endif 

LACK ActionBank 

SACL ISR_TMP 

OUT ISR_TMP,BSR 

OUT NewServo,ACTO 

LAC PwmPeriod,2 

SUB ONE.2 

SUB NewServo 

SACL NewServo 

OUT NewServo,ACTl 

RET 



; !skip end! 
; !END! 



; END OF FILE 



; Set bank for timers 



; Output channel 1 timer 

; Calculate complimentary output 



; Set complimentary channel 2 timer 
; All done., back to the ISR already 
; in progress. 



.end 
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APPENDIX B 
Fuzzy Logic Code 

.title "Fuzzy Servo compensator Revision: 3.4 $" 

SHeader:: C:/src/cl4/ps/vcs/comp.s_v 3.4 01 Oct 1991 17:24:18 "$ 

$config$='YT8 /K! /L;*.ref.def /R:*- /B80" 
! config ! ="/Mcomp.s" 

!NAME! 
compfuz.s 

IPATHS! 

modules 

!0! 

DESCRIPTION 

Servo compensator - uses a fuzzy algorithm (modifying PID code) 

PRINCIPALAUTHORS: 
Dave Sewhuk 

Joe George w/ Fuzzy based on C++ code by Phillip Thrift 
CREATION DATE: 

December 23, 1990 22:54:08 

Jan. 1992 
COPYRIGHT NOTICE: 

(C) Copyright 1990 Teknic Inc. All rights reserved. 

(C) Copyright 1992 Texas Instruments Incorporated. All rights reserved. 

lend! 

*************************************************************************** 

*************************************************************************** 

HEADERS UTILIZED 

.include "macrodef.inc" 
.include "cl4io.inc" 

*************************************************************************** 
*************************************************************************** 

!NAME! 

comp.s 



IPATHS! 
!0! 



Exported Variables 

.def PwmChannel,PwmPeriod 
.def NewServo,ErrDiff,Kintegrator 
.def ErrNow,ErrLast,DesiredPosition 
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.def ninput, noutput, triangle, nmemb 

.def NM, NS, ZE, PS, PM, nloc, NMw, NSw, Zw, PSw, PMw 

.def nrule, rule_table, input, output, rule, memb 

.def loc, num. den, TEMP, TEMPH, x, X, U 



.bss PwmChannel.l 

.bss PwmPeriod,l 

.bss NewServo.l 

.bss ErrDiff,l 

.bss Kintegrator,l 

.bss ErrNow.l 

.bss ErrLast,l 

.bss DesiredPosition,! 



; Current PWM channel 

; Current period of PWM channels 

; New servo value 

; Error difference 

; Error integrator 

; Current error function 

; Last error function result 

; Servo to this position 



; The following are fuzzy section, initialized and uninitialized 



scalein .set 256/3 

scaleout .set 2*100/3 

; .sect "pfuzzin" 

.asect "pfuzzin", lOh 
.label prule_table 



rule table: 



.word 2,2,2 
.word 3,2,1 
.word 4,2,0 
.word 1,2,3 
.word 0,2,4 
.word2,l,3 
.word2,0,4 
.word2,3,l 
.word 
.word 3, 1,2 
.word 1,3,2 



2,4,0 



ninput 
noutput 



.word 2 
.wordl 



triangle .word 3 

nmemb .word 5 

NMaddr .word NM 

NM .word -3*scalein, -2*scalein, -l*scalein 

NS .word -2*scalein, -l*scalein, 0*scalein 

ZE .word -l*scalein, 0*scalein, l*scalein 

PS .word 0*scalein, l*scalein, 2*scalein 

PM .word l*scalein, 2*scalein, 3*scalein 



nloc 

NMwaddr 

NMw 

NSw 

Zw 

PSw 

PMw 



.word 5 
.word NMw 
.word -2*scaleout 
.word -l*scaleout 
.word 0*scaleout 
.word l*scaleout 
.word 2*scaleout 



nrule 



.word 11 



rule tableaddr .word rule table 



input 

output 

rule 

memb 

loc 

num 

den 

minZ 



.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 



TEMP 

TEMPH 

ONE2 



.word 
.word 
.word 1 



xaddr 
Xaddr 
Uaddr 



.word x 
.word X 
.word U 



rule_end: 

rambeg 

x 
X 

u 



.usect "fuzzin", rule_end-rule_table 

.usect "fuzzunin", 2 
.usect "fuzzunin", 2*11 
.usect "fuzzunin", 5 



; lend! 

.**»***»****************»*************************************************** 



*************************************************************************** 



1NAME! 



1PATHS! 



comp.s 



Local Variables 



!0! 



; lend! 



*************************************************************************** 
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.a!*************************************************************************** 



.text 

.def servo ISR 



; !NAME! 
servo_ISR: 

1PATHS! 



functions\all 

functions\interrupt 

modules\utils.s 

!0! 

DESCRIPTION: 

This is the servo compensator. The distribution algorithm is 
modified from PID to fuzzy format. 

RETURNS: 

No changes to foreground environment. 

STACK LEVELS: 


EXTERNAL DATA REFERENCES 

.ref Position ; Position from feedback ISR 

.ref ONE.CONSTFFFF ; Useful constants 

.ref ISR_TMP ; Scratch 

************************************************************** 
iskip start! 

SOVM ; Saturate math 

; First need to copy data over from program 

LACK prulejable 
ADD ONE, 11 
ADD ONE, 8 

LDPK 
LARK ARO, 6Fh 

LARK AR1, rambeg ; beginning of RAM block 

RAMLOOP: MAR *,1 

TBLR *+ 
LDPK 1 
ADD ONE 
LDPK 
MAR*,0 

BANZ RAMLOOP 

; Get Theta and dTheta values 



LDPK 

LAC ErrNow 
SACL 

ZALH 
SUBH 
SACH 

LDPK 

MAR 
LAR ARO.xaddr 
SACH *+ 

LDPK 
LAC ErrNow 
SUB ErrLast 
SACL ErrDiff 

LDPK 
SACL 



ErrLast 

DesiredPosition 

Position 

ErrNow 



*,0 




* 



; Calculate X 

; indices for X[input][memberfn] 
xfn: 

ZAC 
SACL 
xfnil: 

SACL memb ; for (j=0; 

LAC 
SUB 
BLEZ 
xfnml: 
SUB 
BLEZ 
LAC 
ADD 
SACL 
LAR 
LAC 
LT 
MPY 
APAC 
SACL 
LAR 
CALL 



SACL 

LT 

MPY 

PAC 

ADD 



input 
ZAC 
memb 
ninput 
input 
ilend 

LAC nmemb 
memb 
mlend 
xaddr 
input 
TEMP 
ARO, TEMP 
NMaddr 
triangle 
memb 

TEMP 
AR1, TEMP 
membership 

TEMPH 

nmemb 

input 



; Save old error 
; Calc new error 

; store x input value (Theta) 



; dTheta 



; for (i=0; 

; i<ninput 

; j<nmemb 

; load &x[input] ARO 

; load &memb into AR1 



; Get membership value with AR's loaded 
; and leave answer in ACC 
; Save value 

; Calculate index for X[input][memb] 
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ADD 


Xaddr 


SACL 


TEMP 


MAR 


M 


LAR 


AR1, TEMP 


LAC 


TEMPH 


SACL 


* 


LAC 


memb 


ADD 


0NE2 


SACL 


memb 


B 


xfnml 


mlend: 


LAC input 


ADD 


0NE2 


SACL 


input 


B 


xfnil 


ilend: 


NOP 


; exit input loop 


; Zero out U values 


LAR 


ARO, nmemb 


LAR 


AR1, Uaddr 


ZAC 




Uzl:MAR 


M 


SACL 


*+ 


MAR 


*,o 


BANZ 


Uzl 


; calculate value of Z 


zfn: 




ZAC 




SACL 


rule 


zfhrl: 


ZAC 


SACL 


input 


LAC 


nrule 


SUB 


rule 


BLEZ 


zrlend 


; LAC 


ONE2 


LAC 


PM 


SACL 


minZ 


zfnil: 


LAC ninput 


SUB 


input 


BLEZ 


zilend 


LAC 


ninput 


ADD 


noutput 


SACL 


TEMP 


LT 


TEMP 


MPY 


rule 


PAC 




ADD 


input 



; Store value in X[input][memb] 



; memb++ 



; input++ 

; exit memb loop 



; for (i=0; 

; for Q=0; 
; i<nrule 



; initialize Z[rule] 

; kluge since PM = l*scalein 

; i<ninput 



; find &RULE_TABLE[rule] [input] 



ADD 


rule_tableaddr 


SACL 


TEMP 


MAR 


*,o 


LAR 


ARO, TEMP 


LAC 


* 


SACL 


TEMP 


LT 


nmemb 


MPY 


input 


PAC 




ADD 


TEMP 


ADD 


Xaddr 


SACL 


TEMP 


LAR 


ARO, TEMP 


MAR 


\o 


T AC 


minZ 


SUB 


• 


BLZ 


notmin 


LAC 


• 


SACL 


minZ 


notmin: 


LAC input 


ADD 


ONE2 


SACL 


input 


B 


zfnil 


; check max for U 


zilend: LAC 


ninput 


ADD 


noutput 


SACL 


TEMP 


LT 


TEMP 


MPY 


rule 


PAC 




ADD 


ninput 


ADD 


rule_tableaddr 


SACL 


TEMP 


MAR 


*,o 


LAR 


ARO, TEMP 


LAC 


* 


ADD 


Uaddr 


SACL 


TEMP 


MAR 


*,o 


LAR 


ARO, TEMP 


LAC 


minZ 


SUB 


* 


BLZ 


notmax 


LAC 


minZ 


SACL 


* 


notmax: 


LAC rule 


ADD 


ONE2 



; find RULE_TABLE[rule][input] 



; TEMP = RULE_TABLE[rule][input] 
; find &X 



;&(X[input][RULE_TABLE[rule][input]) 
; check for min 



; input++ 



; get &U[RULE_TABLE[rule][input+output] 



; total kluge for input+output 



; got RULE_TABLE[rule][input+output] 



;U[RULE_TABLE[rule][input+output]] 
;if U<minZ (looking for max) 

store new U 

; rule++ 
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SACL rule 
B zfnrl 
zrlend: NOP 



; Need to defuzzify output, u = num/den = sigma U*Loc/ sigma U 
; (ignoring wt.) 
defuzz: 
ZAC 

SACL loc ; for (i=0; 

SACL 
SACL 
LAR 
LAR 



; Sigma loc*U 

Ml: 
SUB 
BLEZ 
MAR 
LT 
MAR 
MPY 
ZALH 
ADDS 
APAC 
SACH 
SACL 
LAC 
ADD 
SACL 
B 

dlendl: LAC 
SACL 
NOP 

; Sigma U 
ZAC 
SACL 
SACL 
LAR 
MAR 

dfl2: 

SUB 

BLZ 

LAC 

MAR 

ADD 

SACL 

LAC 



loc 

TEMP 
TEMPH 
ARO, NMwaddr 
AR1, Uaddr 

LAC nmemb 
loc 

dlendl 

*,0 

*+ 

M 
*+ 

TEMPH 
TEMP 

TEMPH 

TEMP 

loc 

ONE2 

loc 

dfll 

TEMP 

num ; Store 



loc 

TEMP 
ARO, Uaddr 
*,0 

LAC nloc 
loc 

dlend2 
TEMP 
*,0 
*+ 

TEMP 
loc ; ++i 



; i < nloc; 



;i++ 



; for (i=0; 



; i<nloc; 



ADD 


ONE2 


SACL 


loc 


B 


dfl2 


dlend2: LAC 


TEMP 


BZ 


dzero 


SACL 


den 


NOP 




B 


cdiv 


dzero: 


ZAC 


B 


endfuz 



; divide num/den and scale for new servo 



cdiv: 



LT 

MPY 

PAC 

SACH 

LAC 

ABS 

SACL 

ZALH 

ADDS 

ABS 

LARK 



DIV: 



BANZ 

SACL 

LAC 

BGEZ 

ZAC 

SUB 

SACL 

done: 
LAC 
B 

membership: 
MAR 
LAC 
MAR 
SUB 
BLEZ 
MAR 
LAC 
MAR 
SUB 
BLZ 
MAR 



MAR*,0 
TEMPH 
den 

TEMP 
den 

den 

TEMPH 
num 

AR0.15 

SUBC 

DIV 

TEMPH 

TEMP 

done 

TEMPH 
TEMPH 

ZAC 

TEMPH 

endfuz 



*,0 
* 



M 



*+ 

zvalue 

*,0 
* 

M 
*+ 

divl 
*,0 



; Sign 



den 



; Quotient 
;Sign 



; New Servo 



; &x[input] 

; &memb 
; if x>xO 



; if x>xl 
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LAC 


• 




MAR 


M 




SUB 


*+ 


; if x>x2 


BUZ 


div2 




B 


zvalue 




divl: MAR 


*,AR0 


; x-xO 


LAC 


• 




MAR 


*,AR1 




MAR 


*_ 




MAR 


*_ 




SUB 


* 




RET 






div2: MAR 


*,AR1 


; x2-x 


MAR 


*_ 




LAC 


* 




MAR 


*,AR0 




SUB 


* 




RET 






rvalue: LACK 





; value = 


RET 






endfuz: LDPK 


1 




SACL 


NewServo 


; Save result 


ZAC 




; negation kluge 


SUB 


NewServo 




SACL 


NewServo 




; yoda: B 


servoJSR 





Check the answer so that it fits within the limits of the timer's 
period: [O..PwmPeriod*4) 



LACK 
SACL 
LAC 
SUB 
SUB 
BGZ 
LAC 
SUB 
SACL 
B 

pwmisr20: 
LAC 
ADD 
BGEZ 
ZAC 



OFFh 

PwmPeriod 

PwmPeriod,l 

ONE,3 

NewServo 

pwmisr20 

PwmPeriod, 1 

ONE.3 

NewServo 

pwmisrSet 

NewServo 
PwmPeriod, 1 
pwmisrSet 



Check positive: 2*Period-4 
- value 

Oops, pin at biggest 



; Negative OK? 
; Add largest offset 



SUB PwmPeriod,l 
SACL NewServo 



; Make largest negative 



; Set Output PWM to new value 
pwmisrSet: 



LAC 


PwmPeriod,l 


ADD 


NewServo 


SACL 


NewServo 


.if 


ChipVlRl 


.ref 


CONST3 


SUB 


CONST3 


BLEZ 


pwmisrlO 


LAC 


NewServo 


AND 


CONST3 


SUB 


CONST3 


BNZ 


pwmisrlO 


LACK 


4 


ADD 


NewServo 


SACL 


NewServo 


pwmisrlO: 




.endif 




LACK 


ActionBank 


SACL 


ISR TMP 


OUT 


ISR_TMP,BSR 


OUT 


NewServo,ACTO 


LAC 


PwmPeriod,2 


SUB 


ONE.2 


SUB 


NewServo 


SACL 


NewServo 


OUT 


NewServo,ACTl 


RET 





; Iskip end! 
; !END! 



; END OF FILE 
.end 



; Get center period value 
; Add newly computed value 

; PWM bug workaround 



; No fixes needed 
; Lower 2 bits set? 



; OK, No fixes needed 



; Set bank for timers 



; Output channel 1 timer 

; Calculate complimentary output 



; Set complimentary channel 2 timer 
; All done., back to the ISR already 
; in progress. 
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The Programmable 
Fuzzy Logic Array 



Philip Thrift 
Central Research Laboratories 
Texas Instruments Incorporated 



Introduction 



The Programmable Fuzzy Logic Array (PFLA) produces nonlinear multidimensional mappings by encod- 
ing fuzzy rule systems in a programmable array architecture. It can be used as a component of a graphical- 
user interface (GUI) for designing fuzzy controllers, as well as a software blueprint for mapping onto VLSI 
hardware. An advantage of the PFLA over other fuzzy representations, such as the FAM (Fuzzy Associa- 
tive Memory) [2], is the PFLA's ability to easily visualize several inputs and outputs simultaneously. 



Data Flow 



Figure 1 shows the PFLA data flow. Succeeding text describes PFLA in general mapping terms from P 
inputs to Q outputs, but only P = 2 and Q = 1 are shown in Figure 1. Each item number in the text corre- 
sponds to a step of Figure 1. 
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Figure 1. PFLA Data Flow 
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z. ueimea over eacn mpui range [a;, o,j : i = 1,.., r are iuzzy memoersnip iuncuons pi ft, 

where /!; is the number of membership functions defined for input i. Each membership function 
varies between and 1. The cases shown in Figure 1 are trapezoidal membership functions, 
which are defined in Appendix A. Other parametric families of membership functions can be 

substituted. Also, for each fuzzy membership function pi , there is a corresponding label. A typi- 

i 

cal labeling scheme (labels are not shown in Figure 1) for «,• = 7 is: negative large (NL), negative 
medium (NM ), negative small (NS), zero (ZE), positive small (PS), positive medium (PM), pos- 
itive large (PL). Although the fuzzy sets for each input in Figure 1 appear symmetrically spaced, 
these are not necessarily the optimal settings. If trapezoidal memberships functions are used, 
four numbers r / = jy; ^ j j ^ j must be specified for each membership function pi. 

3. Each input is evaluated by each of its fuzzy membership functions to produce a value jt : 

fl = F\ (x.) : i = 1,..., Pj = 1,..., «j (9 ) 

These values appear in the boxes as shown as a thermometer level, f! also refers to the box that 
contains its value. In the example shown, only two boxes for each input have positive evaluation. 

4. Fuzzy rules are encoded in a crossbar pattern. Corresponding to each of the input boxes ji in 

3, above, is a vertical wire dropping down. A horizontal wire crosses those vertical wires and 
also the vertical wires corresponding to output boxes in 6, below. A connection is indicated by 
a •. For each horizontal wire, there is, at most, one connection per input and output variable. Each 
horizontal connection pattern encodes a rule. For example, the connection pattern on the first 
horizontal line encodes the rule: 

If *i is NS and x 2 is ZE, then is PS. ( 10) 
Four rules are shown in Figure 1. R rules can be specified by a table of numbers: 

r\ /* , s\ ,...s° (ID 

r R .•••> ^ . 4 >••••*« (12) 

where r f is in {l,...,n,,NULL}, y is in {l,...,/wy,NULL}. Here, mj is the number of fuzzy sets 
for output j. The kth horizontal wire specifies the rule 

p-l pi -+ G i G s t (i3) 

The NULL indicates that there is no connection for this input/output variable. In Figure 1, the 
connections would be 

3,4,5 
3,5,4 
4,6,3 
4,5,4 

5. Intercepting the horizontal wires between the inputs and outputs are the A ("wedge") boxes. 
A conventional A operator is the numerical minimum of the values (this is used in Figure 1), 
but other operators are possible (for example, product, or any of a set of so-called t-norms [1]). 
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At the klh A box, this is produced: 

** - A (th-fi >$> m 

If fj is NULL, this argument is omitted. 

6. The V ("vee") boxes are computed according to the values of the connections above them. A 
conventional V operator is the numerical maximum of the values (see Figure 1). The values of 
the V boxes are 

g\ = V [h k : si = j] i = 1,...,G j = 1,. (is) 

Other operators [1] can be substituted for this operator (for example, probabilistic sum: x © y 
= x + y-xy , etc.). 

7. Defined over each output range [c,-,4] : i = h—,Q are fuzzy membership functions _ _ _ > Q m t , 

where wi/ is the number of membership functions defined for output i. Also defined for each out- 
put i is a vector /; of locations quantizing the range: \. = m /«;], where qi is the number of 

locations specified for output i. In Figure 1, 17 locations are shown for output 1. 

w>> = & t (1,), i * 1 Q j = 1,..., m { (16) 

Here, Qi is applied by component to get the resulting vector. The ^ are precomputed and stored 
for each V box. Then, this is computed for output i: 

v, = V* (w j A* g . ; ... )M ,™, A* g]) m 

Here, the "wedge" and "vee" operators are not necessarily the ones used in the boxes above. 
In Figure 1, vj is shown as a sequence of vertical bars; V * is the maximum operation, and A * 
is the product. 

8. The final stage is to compute m from v; for each output: 

Vl ■ 1 

where 1 = [1,...,1] is a vector of Is. 

If, corresponding to each output fuzzy set qi , there is a single distinct location p t ^ = G'fV) , anc * A * 
is the product operation, then the above defuzzification procedure reduces to 

= h? • g) ■ i\ + ... + • g :> • (19) 

W\ ■ g) + ... + H? • g^- 

This provides "weighted point mass defuzzification". In this case, only a "weight" and a location /; 
must be stored for each V box. 

Figure 2 shows the state of the system for particular inputs xi,X2. Thermometer levels indicate the values 
in both the fuzzification and the A and V boxes (here, min and max are the operations performed). 
Darkened rule connections indicate values that propagated through the array (since min and max are used, 
on each output line there is, in general, one max winner and one min winner). Figure 2 also shows the layout 
for a four-input, two-output (4-2) system. In a GUI, you can use point-and-clicks to set the rule connections 
and membership function positioning. 
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Figure 2. PFLA 4-2 Layout 
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Summary 

This is the information that provides the setting for the PFLA: 

#Input fuzzy sets: nj np 

#Output fuzzy sets: mi,...,mQ 
[ai,bi], t \ , i = 1,..., P 

Operators: A V 

[c j ,d ] ],l j ,wj wp, j = 1,..., Q A * 



#Inputs: P 
#Outputs: Q 
Input fuzzy sets 

#Rules: R 
#Rule table: 

Output defuzzifiers 



Appendix A 

A trapezoidal membership function on an interval [a,b] is specified by four numbers ti, t2, t3, t4, satisfying 



a s tj s t2 s t3 s t4 s b 



The trapezoidal function is defined by 



TZ(a,t 1> t2 > t 3> t4 > b)(: 





(x-ti)/(t 2 -t!) 
1 

(t 4 -x)/(t 4 -t3) 




for x in [a,ti] 
for x in (ti,t2) 
for x in [t2,t 3 ] 
for x in (t3,t4> 
for x in [t4,b 



Note that if tj = t2 or = t*, this part of the definition is void. 
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